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I. Nomenclature 


b = propellant dependent coefficient 

c = engine exhaust velocity (km/s) 

C 3 = launch energy (km 2 /s 2 ) 

d = propellant dependent coefficient (km/s) 

/ = state differential equation vector 

g = sea-level gravitational acceleration (m/s 2 ) 

h = initial boundary conditions 

I sp = specific impulse (s) 

J = performance index (yrs) 

K t = tankage fraction 

m = spacecraft mass (kg) 

m net = net mass of spacecraft (kg) 

m pp = mass of power and propulsion system (kg) 

m P Top = propellant mass (kg) 

m tank = mass of tank and propellant feed system (kg) 

N — interior-point state constraints 

P = input power (kW) 

r = spacecraft heliocentric position vector (AU) 

rr-’B = planetary flyby radius (DU) 

r P = target planet heliocentric position vector (AU) 

to = launch date 

tj = arrival date 

i t = planetary flyby date 

T = thrust magnitude (N) 

u = thrust direction unit vector 

x = state vector 

a = specific mass of the power and propulsion system (kg/kW) 

/? = propellant mass-flow rate (kg/s) 
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&.V GA = velocity increment from gravity assist (km/s) 
77 = thruster efficiency 

— terminal state constraints 
6 — longitude angle (deg) 


Subscripts 

0 = initial (launch) 

200 = at heliospheric boundary (r = 200 All) 

/ = final (arrival) 

1 = planetary flyby 


II. Introduction 

The first New Millennium interplanetary space mission proposed by the National Aero- 
nautics and Space Administration (NASA) will use electric propulsion (EP) as the pri- 
mary mode of propulsion . 1 For this particular mission, the spacecraft will use solar electric 
propulsion (SEP) and a 30-cm ion engine with xenon as the propellant. In addition, a fast 
Pluto flyby mission using near-term SEP technology has been recently investigated . 2 Other 
EP mission studies include both manned and unmanned lunar and Mars missions using 
nuclear electric propulsion (NEP ). 3-6 The payload advantages of using low-thrust electric 
propulsion in comparison to conventional chemical propulsion have been demonstrated by 
numerous authors for a variety of space missions . 7,8 

A particular deep-space mission of recent interest is the exploration of the heliospheric 
boundary of our solar system. Mewaldt et al . 9 outlined a mission to send a small spacecraft 
to the heliosphere via chemical propulsion. The heliospheric boundary consists of the 
interaction between the solar wind and interstellar plasma and the upper limit of this 
boundary is estimated to be at about 200 AU . 9 Very little is known about the structure of 
the heliosphere and the primary objective of such a mission would be in situ measurements 
of the particles, plasma, and fields at the heliospheric boundary . 9 In addition, a spacecraft 
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enroutc to the heliospheric boundary may be able to encounter Neptune or Pluto and 
therefore enhance the scientific return. 

This paper presents several mission designs for heliospheric boundary exploration using 
spacecraft with low-thrust ion engines as the primary mode of propulsion. The mission 
design goal is to transfer a 200-kg spacecraft to the heliospheric boundary in minimum time. 
The mission design is a combined trajectory and propulsion system optimization problem. 
Trajectory design variables include launch date, launch energy, burn and coast arc switch 
times, thrust steering direction, and planetary flyby conditions. Propulsion system design 
parameters include input power and specific impulse. Both SLP and NLP spacecraft are 
considered and a wide range of launch vehicle options are investigated. Numerical results 
are presented and comparisons with the all-chemical heliospheric missions from Ref. 9 are 
made. 

III. Mission Definition 

As previously mentioned and in accordance with Ref. 9, the objective of the mission 
is to transfer a spacecraft with a net mass of 200 kg to the heliospheric boundary at 200 
AU within a reasonable total mission lifetime. In order to minimize the total trip time, 
the general direction of the trajectory should be toward the “nose” of the heliosphere, 
which is the direction of the velocity vector of the solar system. In the heliocentric-ecliptic 
coordinate frame, this incoming direction of the interstellar plasma corresponds to 254.5 
deg longitude and 7.5 deg latitude. Figure 1 shows a schematic of the orientation of 
the heliosphere. Therefore, the goal of the trajectory and propulsion system optimization 
problem is to obtain the minimum- time transfer to the heliospheric boundary defined above. 


3 




Figure 1: Orientation of the heliospheric boundary. 

Since the mission definition calls for a net mass of 200 kg, the net spacecraft mass m net 
is defined: 


ninet — TTlprop l^ltank T^pp (f) 

The spacecraft’s net mass represents the usable mass for payload plus the basic spacecraft 
structural mass. The initial mass of the spacecraft m 0 upon injection into heliocentric space 
is computed from launch vehicle performance curves with hyperbolic excess energy (C 3 ) as 
the independent variable. Tank mass is the product of the tankage fraction Kt and 

the total EP propellant mass m prop and the mass of the EP power and propulsion system 
m pp is the product of specific mass a and input power P. The vehicle parameters Kt and 
a represent the assumed technology level for EP spacecraft and their respective values are 
presented in the following sections. 
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IV. Trajectory and Propulsion System Optimization 


Solar Electric Propulsion Model 


The SEP system is modeled with thrust magnitude T dependent upon input power P, 
thruster efficiency r/, and specific, impulse I sp 


T = 


2rjP 


( 2 ) 


where c = gl sp is the engine exhaust velocity and g is the Earth’s sea-level gravitational 
acceleration. Propellant mass-flow rate fi is related to thrust as follows: 



( 3 ) 


It is assumed that xenon is utilized as the propellant for the EP system and therefore 
thruster efficiency rj is determined by the relation below using I sp and propellant-dependent 
coefficients derived from theoretical models and experimental data 4 


r] = bc 2 /(c 2 + d 2 ) (4) 

where b = 0.81 and d = 13.5 km/s. Specific impulse is bounded by the limits shown below 
which are typical for EP systems using xenon 4 

2000 < I sp < 7000 s (5) 

It is assumed that I sp is constant over the entire mission which implies a fixed engine 
operating point with no throttling. 

For SEP spacecraft, power P from the solar arrays diminishes with distance from the 
sun. The power ratio relative to available power at 1 AU is modeled by the performance 
curve presented in Fig. 2 for a solar array with silicon cells. The power is nearly 
proportional to the inverse-square distance from the sun and is considered to be zero for 
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Figure 2: Solar array power ratio relative to 1 AIJ. 

distances greater than 3 AU and for distances less than 0.5 AU due to extreme thermal 
conditions. Specific mass a is fixed at 20 kg/kW for this study which represents near-term 
SEP technology 10 and tankage fraction K t is fixed at 10%. 11 

Nuclear Electric Propulsion Model 

The basic system modeling for an NEP spacecraft is essentially the same as the SEP 
spacecraft model with the exception that thruster input power P is constant and does not 
change with distance from the sun. It is assumed that P does not degrade over time. Since 
the NEP system is assumed to utilize ion thrusters and xenon as the propellant, the basic 
relations for thrust, mass-flow rate, efficiency, and I sp range are identical to Eqs. (2-5). 

In order to utilize the current class of launch vehicles, the nuclear reactor must be in 
the kilowatt range instead of the megawatt range. In Ref. 12, a derivative of the Topaz 
2 reactor is proposed for deep-space interplanetary missions. The improved Stirling-cycle 
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Topaz reactor is estimated to have a total mass of 1405 kg and produce 40 kW of power. 12 
Therefore, the specific mass for this proposed nuclear power system is approximately 35 
kg/kW. In Ref. 5, the specific mass for an ion propulsion system (thrusters plus power 
processors) using a nuclear power source is conservatively estimated at 6 kg/kW. I herefore, 
the specific, mass ft for the power and propulsion system for a kW-range NEP system is 
estimated at 41 kg/kW. Since the specific mass typically decreases with increasing reactor 
size (the so-called “economy of scale” for NEP power systems), a lower limit for power 
output of the Topaz-derived system is set at 30 kW. In other words, the estimated o of 41 
kg/kW is assumed to be applicable only for power output levels above 30 kW. The tankage 
fraction K t remains fixed at 10% for the NEP system. 

Problem Statement 

The optimization problem is presented in this section. The problem statement for the 
heliosphere trajectory and propulsion system optimization problem is as follows: 

Find the initial launch date to, the launch energy Co, the durations of the multiple 
powered arcs, the thrust direction unit vector u(t) for the powered arcs, the changes in 
true longitude for the coasting arcs, the planetary flyby conditions, the input power P 0 at 
1 AU, the specific impulse I sp , and the final arrival date tj which minimize: 

J = tf — to (6) 

subject to the equations of motion in an inverse-square gravity field 

x(t) = f{t,x(t),u(t)) (7) 

with the initial boundary conditions 
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x(t 0 ) = h(t 0 ,c 3 ) = x 0 


( 8 ) 


the interior-point state constraints 


N[x(t % ), t,\ = r(<i) — rp(t,i ) 
the terminal state constraints 


( 0 \ 

0 

W 


ii>[x{tj),tf] = r(tj) - r- 2 oo 
and the inequality constraints (Eq. (5)) and 


( 0 \ 

0 

W 


(9) 


( 10 ) 


m net > 200 kg (11) 

In the above problem formulation, the spacecraft s state is denoted by the 7-element 
vector x(t) where the first six elements are the spacecraft’s orbital elements and the seventh 
element is spacecraft mass. Since the goal is to minimize the total trip time, the performance 
index J is the difference in final arrival date and initial departure date. The equations of 
motion (Eq. (7)) are with respect to a sun-centered, inverse-square gravity field and are 
expressed in terms of a modified equinoctial orbital element set 13 in order to allow for both 
elliptical and hyperbolic heliocentric trajectories. The detailed set of dynamical equations 
of motion for a powered spacecraft is presented in Ref. 13. The initial boundary conditions 
(Eq. (8)) are computed by an accurate solar system ephemeris with the Julian date of exit 
from the Earth’s sphere of influence (SOI) and launch energy C 3 as the input variables. 
Equation (9) represents interior-point state constraints which require matching between 
the spacecraft and target planet position vectors for a planetary gravity assist. The gravity 
assist is modeled as an impulsive velocity change and the details are presented in Ref. 
13. The terminal state constraints (Eq. (10)) require that the final spacecraft position 
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vector r(tj) match the position vector of the heliospheric boundary “nose” at 200 AU, 
r200 = [-52.9904, -191.0773, 26. 1052] r AU, which corresponds to a latitude and longitude 
of 7.5 and 254.5 deg, respectively. 

Solution Method 

Since the trajectory and propulsion system optimization problem involves a mix of dis- 
crete control parameters (launch and arrival dates, powered arc durations, planetary flyby 
conditions, input power, I S p) and continuous control functions (thrust direction unit vector 
u(t) during the powered arcs), the problem outlined above is solved using a direct opti- 
mization method. More specifically, the optimal control problem is replaced by a nonlinear 
programming (NLP) problem which in turn is solved by using sequential quadratic pro- 
gramming (SQP), a constrained parameter optimization method. 14 The SQP code used 
here computes the required partial derivatives with first-order forward differences. 15 The 
thrust direction control functions u(t ) are parameterized by linear interpolation through 
a discrete set of control nodes. Therefore, the control nodes (for each powered arc) are 
SQP design variables. Six nodes for each thrust direction component per powered arc were 
deemed sufficient for acceptable trajectory control. Trajectory generation is accomplished 
by numerically integrating the powered equations of motion (Eq. (7)) with a standard 
fourth-order, fixed-step, Runge-Kutta routine. Acceptable numerical accuracy for the tra- 
jectory was achieved with 200 integration steps per powered arc. Coast arcs are computed 
by using analytical methods from two-body orbital mechanics. Planetary flybys for gravity 
assists are treated as impulsive velocity changes and are computed via two-body analytical 
methods by specifying the flyby orbit plane orientation and flyby periapsis radius (see Ref. 
13 for more details). Therefore, these two parameters (per gravity assist) are also included 
as SQP design variables. Since the optimal control problem is solved as a NLP problem, 
the number of powered and coast arcs as well as the number and sequence of planetary 
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flybys must be specified in advance. 

Finally, the terminal state constraints (Eq. (10)) which maintain arrival at the helio- 
spheric boundary “nose” are enforced through three SQP equality constraints. In addition, 
the interior-point state constraints (Eq. (9)) are also enforced through a set of three SQP 
equality constraints per flyby. The inequality constraints on the I sp range as indicated by 
Eq. (5) are enforced through upper and lower box constraints on the appropriate SQP 
design variable and the payload constraint (Eq. (11)) is enforced by an SQP inequality 
constraint. 

V. Numerical Results for SEP Spacecraft 

There exists a very wide range of trajectory options for the minimum-time transfer 
to the heliosphere boundary. For example, the launch vehicle, number of powered and 
coasting arcs, and number and sequence of planetary gravity assists must be specified in 
advance for use of the direct optimization method. For this reason, the minimum-time 
transfer problem is initially solved for the simplest trajectory cases; that is, cases which do 
not involve planetary gravity assists. By solving these “simple” problems, one can obtain 
a feel for other trajectory options such as launch vehicle selection and number of powered 
and coasting arcs. 

Trajectories without Gravity Assists 

Four different existing launch vehicles are investigated: Med- Lite, Delta II 7925, Atlas 
HAS /Star 48B, and Titan IV/Centaur. The last three are the launch vehicles used in Ref. 
9 and are selected so that a direct comparison can be made between the results presented 
here and the trajectory results using conventional chemical propulsion systems. 

The convergence properties of this numerical problem were found to be improved by 
first solving a maximum-energy problem without constraints on the final trajectory direc- 
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tion. For the maximum-energy problem, it is possible to use simple initial guesses for the 
design variables (such as grouping the burns at perihelion and setting the thrust direction 
u{t) along the local horizontal direction). The maximum-energy problems identify trajec- 
tories that achieve solar-system escape conditions and the corresponding solutions are then 
utilized as the initial guesses for the minimum-time heliopause transfer problem. The fi- 
nal trajectory direction constraints are satisfied during the numerical optimization process 
by adjusting the launch date which essentially rotates the entire trajectory in heliocentric 
space. 

A variety of minimum-time missions are obtained for a range of engine sequence schemes 
and the results are shown in Table 1. Initially, the single-burn solution is obtained for each 
launch vehicle followed by the 2-burn solution, 3-burn solution, and so on. All of the optimal 
missions resulted in a net mass at the lower limit of 200 kg. The optimal launch year is not 
presented in Table 1 since the optimal mission opportunities will repeat annually. The best 
trajectory options without gravity assists exhibit a range of trip times from 46.2 years (Med- 
Lite) to 29.4 years (Titan). In Ref. 9, all orbit transfers using chemical propulsion required 
several planetary gravity assists; the fastest transfers ranged from 34.1 years (Delta) to 21.8 
years (Titan). These gravity-assisted transfers using chemical propulsion will be described 
in further detail in the next section. 


Table 1: Minimum-time missions to the heliospheric boundary for SEP spacecraft without 
gravity assists 


Launch 

vehicle 

Number of 
powered arcs 

Launch 

date 

Trip time, 
yrs 

c 3 , 

km 2 /s 2 

m 0 , 

kg 

Po, 

kW 

Isp-) 

S 

Med-Lite 

2 

Nov 22 

47.3 

1.51 

661.5 

9.4 

4818.6 

Med- Lite 

3 

Jan 28 

46.2 

0.60 

673.7 

8.9 

5012.6 

Delta 

1 

Aug 15 

63.7 

2.23 

1272.9 

22.2 

3288.8 

Delta 

2 

Jan 8 

35.8 

2.78 

1259.5 

18.9 

3982.9 

Atlas 

1 

Dec 23 

46.7 

61.26 

957.7 

19.3 

2608.6 

Atlas 

2 

Feb 9 

32.1 

6.12 

2280.0 

34.2 

3564.4 

Titan 

1 

Jan 2 

37.1 

85.43 

1834.8 

41.3 

2437.8 

Titan 

2 

Mar 14 

29.4 

17.24 

6613.2 

98.3 

3208.2 
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As indicated in Table 1, it is not possible to achieve escape conditions with the Med- 
Lite and a single SEP powered arc. Furthermore, the total trip time decreases as the 
number of burn arcs is increased until an upper limit for the number of burns is reached. 
For the Med-Lite launch vehicle, the fastest transfer to the heliosphere requires 46.2 years 
and involves three burns grouped near perihelion. An attempt to obtain a 4-burn solution 
results in a trip time greater than 46.2 years and subsequently the direct method eliminates 
the fourth burn arc during the optimization process. The upper limit on number ol burns 
is two for the three larger launch vehicles as shown in Table 1. It should be noted that 
increasing the number of perihelion burns will decrease the propellant required for a given 
orbit transfer; however, the heliosphere boundary mission is a minimum-time problem and 
not a minimum-fuel problem. I herefore, including additional burn arcs results in additional 
revolutions about the sun which increases the total trip time. 

Table 1 also shows the optimal trajectory and EP system parameters for the orbit 
transfers for each launch vehicle. The EP power and C 3 both increase and I sp decreases as 
the launch capability is increased. Therefore, the thrust magnitude of the SEP system is 
increased in order to minimize the transfer time as larger launch vehicles are employed. All 
of the minimum-time transfers exhibit perihelion burns near a radial distance of about 0.7 
AU so that the solar array performance (as depicted by Fig. 2) is maximized. Furthermore, 
the final burn switching points are at a radial distance of 3 AU which corresponds to the 
maximum radial distance for solar array output. Figure 3 shows the heliocentric trajectory 
and the two perihelion burns for the Delta launch vehicle. In this case, the trajectory 
achieves solar-system escape conditions during the second burn just prior to perihelion 
passage. 

Trajectories with Gravity Assists 

A variety of planetary gravity assist options are explored for the range of launch vehicles. 
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Figure 3: Minimum-time SEP transfer: Delta without gravity assists. 


A good guess for the optimal number of burn arcs for each launch vehicle is determined 
from the solutions in the previous section. 

The strategy employed here is to obtain optimal missions for the simplest gravity assist 
scenario possible and subsequently solve more difficult problems with increasingly complex 
gravity assist maneuvers. Therefore, the initial optimal trajectories involve a single Jupiter 
gravity assist after all SEP powered arcs and escape conditions have been achieved. This 
problem is readily solved by employing an initial guess for the optimization routine based on 
the previous minimum-time orbit transfers without gravity assists. Furthermore, the launch 
date is adjusted until the escape trajectories from the previous solutions pass near Jupiter. 
Two additional trajectory design variables are included in this problem: orientation of the 
planetary flyby plane and periapsis radius of the flyby. The lower limit of the periapsis 
radius is set at 5 Jupiter radii in accordance with gravity assist analyses. 16 

Table 2 presents the optimal trajectory and SEP system parameters for minimum-time 
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transfers with a Jupiter gravity assist for the four launch vehicle options. The gravity- 
assisted trajectories show significant reductions in total mission time over the trajectories 
without gravity assists and these mission-time reductions are more dramatic for the smaller 
launch vehicles. The respective trip-time reductions from using the Jupiter flyby aro' 11.1 
years (Med-Lite), 5.4 years (Delta), 4.2 years (Atlas), and 3.0 years (Titan). As expected, 
the optimal perijove radius for the gravity assist is 5 Jupiter radii for all cases so that 
the maximum AVqa is achieved. The magnitude of the gravity assist is shown in 'Table 2 
and ranges from 14.8 to 17.7 km/s. The 2-burn solutions for the smaller launch vehicles 
(Med-Lite and Delta) exhibit launch dates in late 2001 and the 2-burn solutions for the 
larger launch vehicles (Atlas and Titan) show launch dates in early 2003. The Med-Lite 
and Delta trajectories require an additional year to reach Jupiter’s orbit compared to the 
Atlas and Titan trajectories. 


Table 2: Optimal trajectory and SEP parameters for minimum-time transfers with a Jupiter 
gravity assist 


Launch 

vehicle 

Number of 
powered arcs 

Launch 

date 

Trip time, 
yrs 

c 3 , 

km 2 /s 2 

a, 

kW 

1 sp 1 
S 

km/s 

Med-Lite 

3 

Dec 24 2009 

37.1 

0.98 

9.23 

4705.5 

17.2 

Med-Lite 

2 

Nov 1 2001 

35.1 

1.63 

9.26 

4780.4 

17.7 

Delta 

2 

Dec 20 2001 

30.4 

3.29 

18.84 

3964.5 

16.1 

Atlas 

2 

Jan 17 2003 

27.9 

7.86 

34.82 

3477.7 

15.4 

Titan 

2 

Mar 2 2003 

26.4 

17.92 

96.87 

3220.6 

14.8 


Initially, the minimum-time transfer with the Med-Lite was found with three perihelion 
burns and the result is a trip time of 37.1 years with a launch date in late 2009. The 
launch date for all cases is constrained to be greater than 1999, which is a reasonable time 
frame for SEP technology. The 3-perihelion burn case requires an extra orbit about the 
sun which adds 3.8 years to the time required to reach Jupiter’s orbit as compared to the 
2-burn case. Therefore, in order to flyby Jupiter, the 3-burn trajectory must depart Earth’s 
orbit 3.8 years earlier than the 2-burn case (which is unacceptable since the 2-burn launch 
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date is late 2001) or delay the launch date about 8 years since Jupiter’s period is 11.86 
years. Recall that for the Med-Lite trajectories without gravity assists, the 3-burn case only 
slightly reduced the trip time to the heliospheric boundary compared to the 2-burn case. 
Furthermore, the energy of the ballistic trajectory at Jupiter’s orbit is only slightly greater 
for the 3-burn case as compared to the 2-burn case. Therefore, when a Jupiter gravity assist 
is included, the total mission time is substantially reduced by utilizing a 2-burn trajectory 
which reaches Jupiter’s orbit over 3 years sooner than the 3-burn case. 

Next, a powered Earth gravity assist is included during the final perihelion burn along 
with the coasting Jupiter gravity assist. The lower limit of tin* perigee flyby altitude is set 
at 350 km. Table 3 shows the mission and SEP optimization results for the four launch 
vehicles. All four cases utilized a 2-burn trajectory and all total trip times were reduced 
with the largest reduction at 3.3 years (Med-Lite) and the smallest at 1.2 years (Atlas). 
Figure 4 shows the heliocentric trajectory for the SEP spacecraft using the Delta launch 
vehicle. Again, the trajectory achieves escape conditions during the second burn prior to 
the Earth flyby. The additional gravity assist results in an increase in C 3 and P for all four 
launch vehicles. 


Table 3: Optimal trajectory and SEP parameters for minimum-time transfers with Earth- 


Jupiter gravity assists 1 


Launch 

vehicle 

Launch 

date 

Trip time, 
yrs 

c 3 , 

km 2 /s 2 

Pi b, 
kW 

I Spl 

S 

Total AVga, 
km/s 

*(</), 

deg 

Med-Lite 

Nov 26 2002 

31.8 

2.22 

9.36 

4469.9 

21.5 

281.0 

Delta 

Dec 25 2002 

28.4 

4.21 

19.40 

3708.3 

19.5 

271.7 

Atlas 

Jan 12 2003 

26.7 

9.52 

35.66 

3351.7 

18.6 

267.5 

Titan 

Jan 31 2003 

24.8 

22.43 

98.90 

3103.4 

18.0 

235.3 


1 A11 solutions have two perihelion powered arcs 


Note that the EJ trajectory solutions presented in Table 3 do not satisfy the final 
longitude constraint of 0(tj) = 254.5 deg. For these complex trajectories with multiple 
perihelion burns and two planetary flybys, the final longitude constraint was relaxed since 
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Figure 4: Minimum-time SEP transfer: Delta with EJ gravity assists. 

the strict enforcement of the 0(tj) constraint added 3-6 years to the respective trip times. 
The final longitude angles in Table 3 show deviations from the desired 0(tj) ranging from 
—26.5 to 19.2 deg. It is not clear from Ref. 9 if an exact final longitude angle is necessary (or 
enforced in the subsequent solutions) since formal boundary conditions are not presented. 
Furthermore, Ref. 9 states that “the general direction of the trajectory should be toward the 
nose of the heliosphere, which corresponds to the ram direction of the inflowing interstellar 
gas.” Since the trajectories presented in Table 3 are in the “general” direction of 0(tf) = 
254.5 deg, the final longitude constraint is relaxed in favor of substantially reducing the 
trip time. It is interesting to note that when the Jupiter gravity-assisted trajectories (Table 
2) are re-computed with the final longitude constraint relaxed, the resulting trip times are 
only slightly reduced. In particular, the trip time reduction ranges from a minimum of 0.8 
days (Med-Lite) to a maximum of 95.5 days (Titan). The “worst” final longitude angle for 
these relaxed cases is 0(t/) = 239.5 deg (Titan). 
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For comparison purposes, the best minimum-time transfers using chemical propulsion 
from Ref. 9 are presented in Table 4. These trajectories employ multiple flyby sequences 
including Venus- Earth- Earth- Jupiter gravity assists with a powered Solar flyby (VEEJS) 
and Earth- Jupiter-Solar flyby (EJS). Both single stage (7 jp = 308 s) and 2-stage (I sp = 
‘290 s) chemical propulsion systems are utilized in the analysis in Ref. 9. As indicated by 
comparing Tables 3 and 4, the best SEP trajectory using a Delta exhibits a trip time 5.7 
years shorter than the best chemical propulsion trajectory with the same launch vehicle. 
In addition, the best SEP trajectory with a Med-Lite has a trip time 2.3 years shorter 
than the best chemical propulsion trajectory using a Delta. Furthermore, the trip time 
differences between the SEP and single-stage chemical systems using the Atlas and Titan 
are within one year. 


Table 4: Minimum-time transfers using chemical propulsion (from Ref. 9) 


Launch 

vehicle 

Gravity assist 
sequence 

Trip time (1 -stage), 
yrs 

Trip time (2-stage), 
yrs 

Delta 

VEEJS 

- 

34.1 

Atlas 

VEEJS 

27.7 

25.4 

Titan 

EJS 

24.5 

21.8 


It should be noted that an initial Venus gravity assist and multiple Earth gravity assists 
were attempted with the SEP trajectories since this strategy seems optimal for many of 
the chemical propulsion trajectories. In all flyby cases (VEJ, VEEJ, or EEJ), the SEP 
trajectories showed poorer performance than the EJ trajectories presented in Table 3. A 
potentially attractive gravity assist option is a double Jupiter flyby, this option is not 
investigated here and is reserved for future work. 

Trajectories with Earth-Escape Spirals 

In this section, the hyperbolic Earth-escape phase via the upper stage chemical boost 
from the respective launch vehicle is replaced by an SEP-powered Earth-escape spiral tra- 
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jectory from circular Earth orbit. Since long-term exposure to the Van Allen radiation 
belts can cause substantial performance degradation to the solar arrays, the initial circu- 
lar high Earth orbit (HEO) is above the radiation belts with an altitude of 51,024 km. 
Again, the initial spacecraft mass in HEO is determined by the respective launch vehicle 
performance curves. The SEP spiral time from HEO to Earth -escape conditions (63 = 0) 
is computed by an analytical expression developed by Perkins for universal low-thrust tra- 
jectory solutions . 17 These universal solutions can be scaled by specifying the initial circular 
orbit radius, initial thrust acceleration, and gravitational parameter of the attracting body. 
The minimum-time transfers to the heliosphere (using an EJ gravity assist) for the four 
launch vehicles are presented in Table 5. It can be seen that the trajectories utilizing an 
SEP-powered Earth-escape spiral show very little difference in trip time when compared 
to the trajectories using a chemically-boosted hyperbolic escape trajectory. 1 he lack of 
enhanced performance from the addition of the escape spiral can most likely be attributed 
to the relatively short spiral escape time to C 3 = 0. Since the spiral escape is short (due 
to the high-energy HEO), the advantages of the SEP system are diminished. In the case 
of the Titan, the increase in P (due to including the escape spiral) does not offset the high 
launch energy capability and the total trip time is actually increased by 0.7 years when the 
escape spiral is implemented. 

Table 5: Optimal trajectory and SEP parameters for minimum-time transfers using an 

Earth-escape spiral to C 3 = 0 1 

Launch Launch Trip time, Escape time, P 0 , h P , Total A V G a, #(*/), 


vehicle 

date 

yrs 

days 

kW 

S 

km/s 

deg 

Med-Lite 

Sept 8 2002 

31.3 

51.8 

11.15 

4592.2 

21.4 

251.7 

Delta 

Oct 8 2002 

28.1 

43.6 

22.54 

3945.1 

19.8 

243.1 

Atlas 

Oct 20 2002 

26.7 

41.0 

42.15 

3657.8 

18.9 

239.3 

Titan 

Oct 29 2002 

25.5 

39.9 

140.1 

3394.9 

18.3 

236.3 

'All solutions employ EJ gravity assists and have two perihelion powered arcs 
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Sensitivity Analysis of SEP and Gravity Assist Parameters 

The sensitivity of the performance index ./ (total trip time) to changes in SEP tech- 
nology parameters (o: and I< t ) and a gravity assist parameter (Jupiter flyby radius r FB ) 
was computed by simple finite-difference methods. The goal is to estimate the sensitivity 
OJ/dX where X is the independent parameter of interest (a, /\ t , or r FB ). The respective 
sensitivities were computed by individually perturbing the parameters by 5% from their 
fixed values and re-computing the minimum-time transfers. Recall that the nominal values 
for these parameters are o = 20 kg/kW, K t = 10%, and r FB — 5 ftj ov . The 2-burn, EJ 
gravity-assisted, minimum-time transfers from Table 3 were utilized as the baseline tra- 
jectories for the sensitivity analysis. Since the Jupiter gravity assist provides the largest 
A Vq,\ contribution, only that flyby radius was perturbed for the analysis. 

Table 6 presents the three performance sensitivities for the four launch vehicles. In gen- 
eral, the performance of the smaller launch vehicles (Med- Lite and Delta) is more sensitive 
to changes in o and t fb and less sensitive to changes in than the perfoimance of the 
larger vehicles. However, the performance sensitivity due to the SEP parameters a and A t 
exhibits a slight variation across the launch vehicle spectrum. In contrast, the performance 
sensitivity due to the Jupiter flyby radius nearly doubles between the Med-Lite and Titan 
vehicles. Using this simple linear analysis, increasing a from 20 to 35 kg/kW would (on 
average) increase the trip time by 4.9 years and increasing K t from 10% to 20% would only 
increase the trip time by one year. Furthermore, increasing the minimum r FB from 5 to 
10 Jupiter radii would increase the trip time by a range of 1.8 years (Titan) to 3.5 years 
(Med-Lite). 

VI. Numerical Results for NEP Spacecraft 

The minimum-time transfers to the heliospheric boundary were obtained for the NEP 
spacecraft outlined previously and the results are presented in Table 7. Trajectories utilizing 
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Table 6: Performance sensitivity to SEP and gravity assist parameters 1 


Launch 

vehicle 

dJl da, 
yrs-kW/kg 

dJ/dK u 

yrs/% 

dJ /dr fb, 

yrs/ R Jov 

Med-Lite 

0.36 

0.09 

0.70 

Delta 

0.32 

0.10 

0.49 

Atlas 

0.31 

0.10 

0.41 

Titan 

0.31 

0.12 

0.36 


1 All trajectories employ EJ gravity assists 


a chemically-injected Earth escape and an NEP-powered spiral Earth escape are presented 
in this table. The initial circular Earth orbit for the spiral Earth escape is set at a ‘‘nuclear 
safe” altitude of 800 km. 3 All solutions employed a single Jupiter gravity assist and a 3-burn 
engine sequence was assumed for all cases. 


Table 7: Optimal trajectory and NEP parameters for minimum-time transfers with a 


J upiter gravity assist 


Launch 

vehicle 

Launch 

date 

Trip time, 
yrs 

Escape time, 
yrs 

C3, 

km 2 /s 2 

P, 

kW 

I spi 
S 

AVfcvi, 

km/s 

Delta 1 

Mar 6 2000 

22.1 

1.5 

0.0 

30.00 

7000 

18.8 

Atlas 

Oct 1 2001 

25.8 

0.0 

0.63 

30.00 

7000 

18.8 

Atlas 1 

Mar 18 2000 

21.5 

1.4 

0.0 

50.14 

7000 

15.6 

Titan 

Dec 22 2000 

18.2 

0.0 

22.81 

37.56 

7000 

18.7 

Titan 1 

Apr 1 2000 

20.9 

1.4 

0.0 

152.32 

7000 

16.0 


1 Utilizes Earth-escape spiral to 63 = 0 


The resulting NEP trajectories are much more “direct” than the SEP trajectories since 
NEP power is constant and the need to perform perihelion burns is eliminated. Therefore, 
the long aphelion coasting arcs are either removed altogether during the optimization pro- 
cess or replaced by short coasting arcs (Ad ~ 30 deg). These “direct” NEP trajectories 
typically involve near-continuous thrusting out to heliocentric distances ranging from 35 
to 40 AU. As an example of these “direct” trajectories, the optimal heliocentric NEP tra- 
jectory for the Delta is presented in Fig. 5. As shown in Fig. 5, the spacecraft essentially 
spirals away from the sun and eventually achieves solar-system escape conditions just prior 
to the Jupiter flyby. The NEP-powered burn arc for this case is terminated at a distance 
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Figure 5: Minimum-time NEP transfer: Delta with Jupiter gravity assist. 

of 36.3 AU (which is not shown in Fig. 5) at which point the spacecraft’s speed is 12.5 
AU/yr. 

As indicated in Table 7, an NEP trajectory could not be obtained for the Med-Lile 
since the smallest possible reactor and engine mass (m pp = 1230 kg) is about 60% of the 
mass delivered by the Med-Lite to the nuclear safe circular Earth orbit. In addition, an 
NEP trajectory using a direct injection into heliocentric space via the upper stage of a. 
Delta could not be obtained since the minimum m pp is about 90% of the injected mass 
with C 3 = 0. The only possible NEP trajectory off a Delta must employ the NEP-powered 
Earth-escape spiral. The Earth-escape spiral lasts about 1.5 years for all three launch 
vehicles and it should be noted that including the escape spiral reduces the trip time for 
the Atlas but increases the trip time for the Titan. The best NEP trajectories for the three 
launch vehicles exhibit better performance compared to the respective SEP cases with trip 
times ranging from 22.1 (Delta) to 18.2 years (Titan). Although an NEP trajectory off 
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a Titan can reach the heliospheric boundary in less than 20 years, the performance of an 
NEP spacecraft off a Delta is comparable and therefore probably preferable. 

VII. Conclusions 

A new deep-space mission to the heliospheric boundary using ion propulsion spacecraft 
has been analyzed. The mission design is accomplished by solving a complex trajectory 
and propulsion system optimization problem. The goal is to minimize the total trip time 
to the heliospheric boundary; the design variables include launch date, launch energy, burn 
times, thrust vector direction, planetary flyby conditions, engine power, specific impulse, 
and arrival date. 

Several optimal mission designs are obtained for spacecraft propelled by solar electric 
propulsion (SEP) and nuclear electric propulsion (NEP) for a wide range of launch vehicles. 
The best SEP trajectories exhibit total trip times ranging from 31.3 years (Med- Lite launch 
vehicle) to 24.8 years (Titan). The best NEP trajectories show trip times ranging from 
22.1 years (Delta) to 18.2 years (Titan). The trip times for electric propulsion spacecraft 
are comparable or less than the trip times for spacecraft propelled by chemical systems. In 
particular, this analysis shows that SEP is capable of performing the mission off smaller 
launch vehicles (Med-Lite) in less time than an all-chemical system off larger launch ve- 
hicles (Delta). The NEP spacecraft show significant reductions in trip time compared to 
an all-chemical system for all launch vehicles. In addition, the overall mission analysis 
shows that the performance gap between large and small launch vehicles is narrowed by 
utilizing either SEP or NEP spacecraft. These results are perhaps surprising since electric 
propulsion has long been identified as an efficient means for performing space missions in 
terms of payload capabilities but not necessarily in terms of trip time savings. This analysis 
also demonstrates the advantages of using low-thrust spacecraft for high-energy deep-space 
missions. 
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